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'  1 1976) .  Configuration  geometry,  attitude,  and  Mach  range  capabilities  are  con¬ 
sistent  with  those  accommodated  by  the  Datcom.  The  program  contains  a  trim 
option  that  computes  control  deflections  and  aerodynamic  increments  for  vehicle 
trim  at  subsonic  Mach  numbers,  -4/olume  I  is  the  user’ 8  manual  and  presents 
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program  capabilities,  input  and  output  characteristics,  and  example  problems. 
Volume  II  described  program  implementation  of  program  Implementation  of  Datcom 
methods^  Volume  III  discusses  a  separate  plot  module  for  Digital  Datcom. 

— , 

The  program  is  written  in  ANSI  Fortran  IV.  The  primary  deviations  from 
standard  Fortran  are  Namelist  input  and  certain  statements  required  by  the  CDC 
compilers.  Core  requirements  have  been  minimized  by  data  packing  and  the  use 
of  overlays. 

User  oriented  features  of  the  program  include  minimized  input  requirements, 
input  error  analysis,  and  various  options  for  application  flexibility. 

/ 


FOREWORD 


This  report,  "The  USAF  Stability  and  Control  Digital  Datcomi"  describes 
the  computer  program  that  calculates  static  stability,  high  lift  and  control, 
and  dynamic  derivative  characteristics  using  the  methods  contained  in  Sec¬ 
tions  4  through  7  of  the  USAF  Stability  and  Control  Datcom  (revised  April 
1976).  The  report  consists  of  the  following  three  volumes: 
o  Volume  I,  Users  Manual 

o  Volume  II,  Implementation  of  Datcom  Methods 
o  Volume  III,  Plot  Module 

A  complete  listing  of  the  program  is  provided  as  a  microfiche  supplement. 

This  work  was  performed  by  the  McDonnell  Douglas  Astronautics  Company, 
Box  516,  St.  Louis,  Missouri  63166,  under  contract  number  F33615-77-C-3073 
with  the  United  States  Air  Force  Systems  Command,  Wright-Patterson  Air  Force 
Base,  Ohio.  The  subject  contract  was  initiated  under  Air  Force  Flight 
Dynamics  Laboratory  Project  8219,  Task  82190115  on  15  August  1977  and  was 
effectively  concluded  in  November  1978.  This  report  supersedes  AFFDL 
TR-73-23  produced  under  contract  F33615-73-C-3058  which  extended  the  program 
to  include  Datcom  Sections  6  and  7  and  a  trim  option;  and  AFFDL-TR-76-45  that 
incorporated  revisions  and  user  flexibility  under  contract  F33615-75-C-3043. 
The  recent  activity  generated  a  plot  module,  updated  methods  to  incorporated 
the  1976  Datcom  revisions,  and  provide  additional  user  oriented  features. 
These  contracts,  in  total,  reflect  a  systematic  approach  to  Datcom  automation 
which  commenced  in  February  1972.  Mr.  J.  E.  Jenkins,  AFFDL/FGC,  was  the  Air 
Force  Project  Engineer  for  the  previous  three  contracts  and  Mr.  B.  F.  Niehaus 
acted  in  this  capacity  for  the  current  contract.  The  authors  wish  to  thank 
Mr.  Niehaus  for  his  assistance,  particularly  in  the  areas  of  computer  program 
formulation,  implementation,  and  verification. 

Requests  for  copies  of  the  computer  program  should  be  directed  to  the 
Air  force  Flight  Dynamics  Laboratory  (FGC).  Copies  of  this  report;  can  be 
obtained  from  the  National  Technical  Information  Serclce  (NTIS) . 

This  report  was  submitted  in  April  1979. 
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SECTION  1 


INTRODUCTION 


In  preliminary  design  operations,  rapid  and  economical  evaluation  of 
aerodynamic  stability  and  control  characteristics  are  usually  essential.  The 
USAF  Stability  and  Control  Digital  Datcom  provides  estimates  for  these  char¬ 
acteristics;  however,  graphical  representation  of  data  is  normally  required 
for  an  expedient  evaluation  of  the  vehicle  aerodynamic  characteristics. 

This  report  describes  the  USAF  Digital  Datcom  Plot  Module.  This  module 
provides  the  user  with  a  graphical  representation  of  data  calculated  by 
Digital  Datcom.  The  Plot  Module  was  developed  as  a  separate  program,  or 
post  processor,  to  provide  the  user  with  maximum  flexibility,  simplicity,  and 
economy.  These  benefits  are  realized  because  the  module  is  not  required  to 
operate  within  the  Digital  Datcom  logic  and  structure,  but  interfaces  with  a 
file  written  from  the  Digital  Datcom  Output  Modules.  A  schematic  of  the  Plot 
Module  interfaces  is  presented  in  Figure  1. 

The  Digital  Datcom  Plot  Module  produces  plots  of  the  following  Digital 
Datcom  data: 

1.  static  stability  characteristics,  including  the  longitudinal 
and  lateral-directional  stability  derivatives; 

2.  wing  downwash  and  dynamic  pressure  ratio  data  in  the  region  of  the 
horizontal  tail,  for  configurations  with  a  wing  and  horizontal 
tail ; 

3.  high-lift  and  control  characteristics  generated  from  symmetrical  and 
asymmetrical  flap  analysis;  ar.d 

4.  trim  data  for  configurations  with  flap  or  all  moveable  horizontal 
tall  controls. 

The  Digital  Datcom  Plot  Module  was  designed  to  operate  in  a  batch  mode 
utilizing  an  "off  line"  plotter.  An  "off  line”  plotting  system  is  one  in 
which  the  program  writes  plot  data  and  commands  to  an  external  file,  usually 
a  magnetic  tape.  This  tape  is  then  used  to  drive  the  plotter  as  Illustrated 
in  Figure  1.  An  "on  line"  version  of  the  program  may  be  generated  by  con¬ 
verting  the  CALCOMP  software  calls  to  the  appropriate  calls  for  the  system  in 
use.  In  an  "on  line"  system,  the  plot  data  and  commands  are  transmitted 


directly  to  the  plotter,  usually  a  Cathode  Ray  Tube  (CRT).  The  CALCOMP 
subroutines  used  by  the  plot  module  are  discussed  in  the  Appendix. 

This  program  is  written  in  the  FORTRAN  IV  language  for  the  CDC  CYBER  74 
and  175  computers.  The  deviations  from  ANSI  standard  FORTRAN  were  transfer 
on  end  of  file  and  those  statements  required  by  the  CDC  compilers.  The 
program  utilizes  the  CALCOMP  plotter  software  available  at  the  AFFDL,  ASD 
computer  center. 

Direct  all  program  inquiries  to  AFFDL  FGC,  Wright-Patterson  Air  Force 
Base,  Ohio,  45433.  Phone  (513)  255-4315. 


SECTION  2 


PROGRAM  CAPABILITIES 

The  Digital  Datcom  plot  module  displays  aerodynamic  data  calculated 
using  Digital  Daicom.  Plots  can  be  obtained  of  the  static  stability, 
high~l.il  and  control,  and  trim  data.  The  user  is  refered  to  Datcom  and 
Volumes  I  and  II  of  this  report  for  a  discussions  of  the  available  data. 

2.1  PROGRAM  DESCRIPTION 

Digital  Datcom  computations  consist  of  one  or  more  case.',  where  each 
case  consists  of  the  aerodynamic  analysts  of  a  configuration  over  a  range  of 
flight  conditions.  Plot  module  output  also  consists  of  one  or  more  caaes, 
with  each  plot  case  producing  plots  of  a  single  Digital  Datcom  case.  The 
aerodynamic  data  plotted  jy  this  program  are  located  on  a  file  written  by 
Digital  Datcom.  The  structure  of  this  file  is  discussed  in  Section  5.  Data 
from  sources  other  than  Digital  Datcom  may  be  plotted  by  creating  a  file  with 
the  proper  format. 

The  plot  module  can  accept  data  from  one  or  two  aerodynamic  data  files 
in  a  single  case  and  up  to  four  different  files  during  each  program  exerj- 
tlon.  This  permits  the  user  to  generate  plots  snowing  comparisons  between 
two  Digital  Datcom  runs  or  between  Digital  Datcom  output  and  test  data.  If 
two  files  are  read  simultaneously,  they  are  assumed  to  have  the  same  flight 
condition  a"d  re^^ence  data.  The  user  has  the  option  to  select  one  of  three 
basic  codes  for  displaying  each  file.  The  three  codes  a"e:  a  solid  line 
between  data  points,  a  solid  line  with  symbols  at  the  data  points,  and 
symbols  at  the  data  points  and  no  line. 

A  Digital  Datcom  aeiodynamic  data  file  contains  the  static  stability 
portion  of  the  Ideal  Output  Matrix  (10M),  see  Appendix  C  of  Volume  1.  There¬ 
fore,  all  of  the  configuration  build-up  data  calculated  by  Digital  Datcom  are 
available  to  the  plot  module,  the  user  has  the  option  to  generate  configura¬ 
tion  build-up  plots  and  specify  the  component  combinations  that  will  be 
plotted. 

2.2  PLOT  GROUPS 

The  plot  module  can  produce  4/  different  plots  of  Digital  Datcom 
data.  In  ordsr  to  more  efficiently  process  these  plots,  they  have  been 


divided  into  seven  groups  of  related  data.  The  types  of  plots  within  each 
group  a.  e  defined  below. 

o  GROUP  1  -  Static  Stability  as  a  Function  of  Angle  of  Attack, 
o  GROUP  2  -  Drag  Polars  and  Stability  Plots 
c  GROUP  3  -  Downwash  Data 
o  GROUP  A  -  Mach  Effects 
o  GROUP  5  -  Symmetrical  Flap  Data 
o  GROUP  6  -  Asymmetrical  Flap  Data 
o  GROUP  7  -  Trim  Data 

Data  from  plot  groups  1  through  A  will  be  plotted  unless  the  user  specifies 
that  the  group  is  inactive.  Within  each  of  these  four  groups,  the  plots  are 
initially  defined  as  active  or  inactive;  however,  the  user  can  override  the 
preset  values.  Active  plots  are  produced  whenever  the  plot  group  is  active. 
Plot  groups  5  through  7  are  nominally  inactive  and  will  be  plotted  only  at 
the  request  of  the  user.  The  plots  in  these  three  groups  are  nominally 
Inactive  and  must  be  specifically  requested  by  the  user.  A  description  of 
the  individual  plots  in  each  group  and  their  default  activity  status  is  given 
in  Section  3.8. 

2.3  PLOT  FORMAT 

The  basic  plot  format  is  shown  in  Figure  2.  The  first  two  lines  in  the 
label  at  the  top  of  the  plot  are  user  defined.  The  next  three  lines  defining 
the  flight  condition  and  reference  data  are  fixed.  On  Mach  number  plots 
(Group  A),  the  Mach  number  and  Reynolds  number  information  is  not  drawn  and 
the  altitude  is  given  instead.  The  axis  labels  are  defined  in  the  program, 
however,  the  user  has  the  option  to  change  them.  Axis  scaling  is  determined 
internally  by  the  program,  with  a  user  option  available  to  override  the 
internal  scaling  and  define  scale  limits.  The  cot. figuration  build-up  combi¬ 
nations  shown  in  Flgu-e  2,  demonstrate  the  labeling  technique  used  when 
multiple  lines  appear  on  a  plot.  Configuration  build-up  combinations  plotted 
are  also  a  user  option. 


FIGURE  7  SAMPLE  PLOT 


SECTION  3 


INPUT  DATA  DEFINITION 

There  are  two  sets  of  input  data  required  by  the  Plot  Module;  Digital 
Datcom  plot  data,  and  plot  module  control  data.  The  Digital  Datcoa  plot  data 
are  located  on  a  file  generated  by  Digital  Datcoa  when  the  "PL0T"  control 
card  is  present  in  a  case  (see  Volume  I,  Section  3.5).  A  maximum  of  four  of 
these  files  may  be  used  during  any  single  execution  of  the  plot  module  and 
are  read  from  1/4  units  11  through  14.  The  plot  module  control  data  define 
which  plots  will  be  generated  and  their  format.  These  control  cards  arc  read 
from  the  normal  input  stream. 

There  are  15  control  cards  recognized  by  the  plot  module.  These  car<»s , 
like  the  Digital  Datcom  control  cards,  have  the  following  characteristics. 

1.  Each  card  consists  of  &  control  word,  starting  in  column  1,  that 
defines  the  program  function  under  control. 

2.  Control  cards  may  appear  in  any  order  in  the  input  stream,  with  the 
exception  of  the  case  termination  cards.  However,  there  are  some 
limitations  with  regard  to  the  PL0T,  SCALE,  and  LABEL  cards. 

3.  The  control  cards  may  appear  more  than  one  time  within  a  case,  the 
last  appearance  will  be  used. 

4.  If  a  given  control  card  is  not  input,  the  program  defaults  are  used. 
Each  control  card  and  the  corresponding  program  defaults  are  defined 
in  the  following  subsections. 

The  plot  module  control  cards  also  contain  input  data  following  the  control 
word.  These  data  are  formatted  inputs,  and  must  be  located  in  the  correct 
columns  in  order  to  be  properly  read.  Five  additional  rules  for  the  plot 
module  control  card  Inputs  are  given  below. 

1.  Input  data  for  an  integer  variable  does  not  leave  a  decimal  point  r*nd 
is  right  Justified  in  the  card  column  field.  Integer  variables 
start  with  the  letters  I  through  N. 

2.  Input  data  for  real  variables  must  have  a  decimal  point  and  may  be 
located  anywhere  in  the  card  column  field.  Real  variables  start 
with  the  letters  A  through  H  and  0  through  Z. 

3.  When  a  control  card  is  used,  all  of  the  required  input  data  must  be 
defined. 


4.  All  input  variables  are  initialized  prior  to  reading  case  data  and 
are  not  saved  for  the  following  case  unless  requested  by  the  user. 

5.  The  control  word  on  each  card  must  be  followed  by  at  least  one 
blank. 

The  plot  control  cards  are  summarized  in  figure  3  and  discussed  iu  detail 
below.  It  is  suggested  that  the  user  study  the  example  problem  in  Section  4. 

All  input  variables  and  common  arrays  are  initialized  "UNUSED”  (1.0  E-bO 
on  CDC  systems)  prior  to  case  execution.  The,  program  uses  the  value  of 
"UNUSED"  as  a  program  switch  to  determine  if  particular  real  variables  have 
been  defined.  In  the  Plot  Module,  these  tests  are  performed  on  the  user 
defined  scale  limits  (see  SCALE  Card),  and  on  data  in  the  Digital  Datccm  plot 
file. 

To  aid  the  user  in  interpreting  the  input  data  definitions,  the  charac¬ 
ter  "0“  is  the  number  zero  and  the  character  "0"  is  the  fifteenth  letter  in 
the  alphabet. 

3.1  FILE  CARD 

The  FILE  card  specifies  the  files  containing  Digital  Datcom  plot  data 
and  defines  the  case  that  will  be  plotted.  Two  separate  files  may  be  defined 
for  each  plot  case  thereby  providing  plotted  comparisons  between  two  sets  of 
data.  These  two  files  are  designated  the  "A”  and  "B"  files.  File  A  is  the 
primary  file  and  is  used  to  read  the  flight  condition  and  reference  data 
and  set  the  plot  scales.  Care  must  be  exercised  when  two  files  are  being 
plotted,  the  scale  limits  are  defined  using  file  A  and  if  come  data  from  file 
B  falls  outside  these  limits  it  will  not  be  plotted.  File  "A"  must  be 
defined  for  every  plot  case,  file  "B"  is  optional.  The  program  default 
defines  file  "A"  on  1/0  unit  1.1  and  file  "B"  as  not  present.  File  "B"  is 
assumed  to  be  absent  whenever  the  1/0  unit  number  is  outside  the  range  of  11 
to  14.  The  user  may  therefore  "turn  off"  file  "B"  by  setting  the  unit  number 
to  zero.  The  format  of  the  FILE  card  and  definition  of  its  input  data  are 
given  below. 

Column  Code  Explanation 

2-4  FILE  Control  word.  The  four  letters  FILE. 

7  ID  File  identification.  Input  either  the 

letter  A  or  B. 
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V  — 


Column 


Explanation 


11-13 


8-10  IUNIT  File  1/0  unit  number,  right  justified  inte¬ 
ger  in  the  range  11  to  14  (11  <  IUNI  <  14). 

If  the  unit  number  is  outside  this  range  for 
file  B,  it  is  not  read  or  plotted.  Defaults: 

File  A,  IUNIT  -  11; 

File  B,  IUNIT  -  0. 

11-13  ICASE  Digital  DATCOM  case  number  to  be  plotted, 
right  justified  integer  (ICASE  >_  1).. 

Defaults: 

File  A,  ICASE  -  1; 

File  B,  ICASE  -  1. 

3.2  LINE  CARD 

The  LINE  card  controls  the  method  of  drawing  lines  between  Digital 
Datcow.  data  points  and  the  drawing  of  symbols  around  the  data  points.  The 
plot  module  provides  three  options  for  drawing  lines  between  points;  straight 
line  between  points,  a  spllned  curve  between  points,  and  no  line  drawn.  The 
user  also  has  the  option  to  have  symbols  drawn  around  the  data  points  in 
conjunction  with  any  of  the  three  line  codes.  Separate  line  code  options  may 
be  assigned  to  files  A  and  B.  The  format  of  the  LINE  card  and  definition  of 
its  input  data  are  given  below. 

Column  Code  Explanation 

1-4  LINE  Control  word.  The  four  letters  LINE. 

7  ID  File  identification.  Input  either  the  letter 

Aj  or  B. 

8-10  ILINE  Line  code  for  the  designated  file,  right 

justified  integer.  The  three  codes  are: 

ILINE  "  0,  straight  line  between 
points 

ILINE  >,  1 ,  spline  curve 

ILINE  <,  -I,  no  line  drawn. 

Defaults: 

File  A,  ILINE  -  0 

File  B,  ILINE  -  0 


Column  Code  Explanation 


11-13  ISYM  Symbol  code  for  the  designated  file,  right 
Justified  integer.  The  codes  are: 

ISYM  ■  0,  no  symbols  are  drawn 
ISYM  -  j,  ( j  ^  l)t  symbols  are  drawn  at 
every  jth  point. 

Defaults: 

File  A,  ISYM  -  0 
File  B,  ISYM  -  0 

3.3  TITLE  CARDS 

Two  control  cards  are  used  to  specify  the  user  defined  titles  that 
are  drawn  at  the  top  of  the  plot  (see  Figure  2).  The  control  card  TITLE1 
specifies  the  title  drawn  on  the  top  line,  and  the  control  card  TITLE2 
specifies  the  title  drawn  on  the  second  line.  The  user  specifying  only  one 
title  line  may  prefer  to  use  only  the  TITLE2  control  card.  The  format  of  the 
TITLE  cards  and  definition  of  their  input  data  are  given  below. 

Column  Code  Explanation 

1-6  TITLE 1  Control  word.  The  five  letters  and  number 

or  TITLE 1  or  TITLE2. 

TITLE 2 

8^9  NC  The  number  of  characters  in  the  title,  right 

justified  integer  <0  <  NC  <  40).  If  NC  -  0, 
the  title  is  not  drawn.  The  maximum  number 
of  characters  is  40,  counting  inbedded 
blanks.  Default: 

TITLE  1,  NC  -  0; 

TITLE 2,  NC  -  0. 

11-50  title  Text  of  the  title.  The  characters  must  be 
left  justified  in  the  field. 

3.4  BUILD  CARD 

The  BUILD  card  permits  the  user  to  select  the  configuration  components 
and  component  combinations  thdt  will  be  drawn  on  each  plot.  The  control 
word,  BUILD,  must  be  followed  by  at  least  cne  blank  and  the  configuration 
codes.  There  are  four  basic  forms  for  the  configuration  codes.  Regardless 
of  the  form  specified,  data  for  the  complete  configuration  will  always  be 


1-5 

7-80 


BUILD 

Code 


ALL 


2. 


N0N2 


3. 


plotted.  The  format  of  the  BUILD  card  and  definition  of  the  four  forms  for 
the  configuration  codes  are  given  below. 

Column  Code  Explanation ■ 

Control  word.  The  five  letters  BUILD. 

1.  One  or  more  of  the  IOM  configuration 
codes  defined  below.  Each  code  must  be 
separated  by  a  blank  or  a  comma. 

All,  the  three  letters  ALL.  All  of  the 
components  and  component  combinations 
defined  in  the  Digital  DATCOM  run  will 
be  plotted. 

None,  the  four  letters  N0NE.  Only  data 
for  the  complete  configuration  will  be 
plotted. 

No  codes  given  columns  7  to  80  are 
blank.  The  IOM  configuration  codes 
B,  BW,  BWH,  and  BWHV  are  selected. 

Configuration  data  not  calculated  in 
Digital  DATCOM  will  not  be  plotted. 

Default: 

Code  -  NONE. 

The  Ideal  Output  Matrix  (IOM)  consists  of  eleven  configuration  components  and 
component  combinations  that  may  be  plotted.  Each  10M  configuration  is 
assigned  a  one  to  four  character  code  that  is  used  on  the  first  form  of  the 
BUILD  card.  These  codes  are  defined  below. 

IOM  Config.  Code  Definition 


blank  4 . 


B 

W 

H 

VT 

VP 

BW 

BH 

BV 


Body  alone 
Wing  alone 

Horizontal  tail  alone 

Vertical  tall  alone 

Ventral  fin  alone 

Body  plus  Wing 

Body  plus  Horizontal  tail 

Body,  Vertical  tail  and/or  Ventral 

fin 
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BWH  Body.  Wing,  Horizontal  tail 

BWV  Body,  Wing,  Vertical  tail  and/or  Ventral 

fin 

BWHV  Body,  Wing,  Horizontal  tail.  Vertical  tail 

and/or  Ventral  fin 

3.5  MACH  CARD 

The  MACH  card  gives  the  user  the  option  to  specify  the  Mach  numbers 
selected  for  plots  at  constant  Mach  number.  Plots  with  Mach  number  as  the 
independent  variable  are  not  affected  by  this  card.  The  Mach  numbers  to  be 
plotted  are  selected  by  specifying  the  Mach  number  index  in  the  Digital 
Datcom  Mach  list.  The  Mach  list  is  defined  in  Volume  I  of  this  report  under 
Namelist  FLTC0N.  If  the  MACH  card  is  not  input,  plots  will  be  generated  at 
all  the  Mach  numbers  in  the  Digital  Datcom  output.  The  format  of  the  MACH 
card  and  definition  of  its  input  data  are  given  below 


Column 

Code 

Explanation 

1-4 

MACH 

Control  word,  the  four  letters  MACH. 

6-8 

IM(  1) 

9-11 

1M(2) 

12-14 

IM(3) 

15-17 

1M(4) 

18-20 

IM<5) 

21-23 

IM(6) 

Index  to  Mach  number  in  Mach  list. 

24-26 

IM(7) 

right  justified  Integers.  Indices 

27-29 

IM(8) 

may  be  in  any  order.  For  example. 

30-32 

IM<9) 

to  specify  the  first  and  third  Mach 

33-35 

IM(10) 

numbers  the  Mach  card  may  be  coded 

36-38 

IM<11) 

in  any  of  the  following  ways: 

39-41 

1M(12) 

MACH  1  3, 

42-44 

IM( 13) 

MACH  3  1, 

54-47 

IM( 14) 

etc. 

48-50 

IM(15) 

51-53 

IM( 16) 

54-56 

IM( 17) 

57-59 

IM(  18) 

60-62 

IM( 19) 

63-65 

IMC2C) 

The  ALT  card  gives  the  user  the  option  to  specify  the  altitudes  at  which 
Digital  Datcom  data  will  be  plotted.  This  card  is  utilized  with  data  that 
were  computed  with  altitude  looping  options  2  or  3.  These  options  are 
defined  in  Volume  1  (namelist  FLTC0N,  variable  name  L00P).  The  altitudes  to 
be  plotted  are  selected  by  specifying  the  altitude  index  in  the  Digital 
Datcom  altitude  list  (see  Volume  I.  namelist  FLTC0N,  variable  ALT).  The 
format  of  the  ALT  card  and  definition  of  its  input  data  are  given  below. 
Column  Code  Explanation 

1-3  ALT  Control  word,  the  three  letters  ALT. 

6-8  IA( 1) 

9-11  IA(2) 

12-14  IA(3) 

15-17  IA(4) 

18-20  IA(5) 

21-23  IA(6) 

24-26  IA(7) 

27-29  IA(8) 

30-32  IA(9)  Index  to  altitude  in  altitude  list,  right 

33-35  IA(10)  justified  integers.  Indices  may  be  in 

36-38  IA(ll)  any  order.  See  the  example  given  with 

39-41  IA(12)  the  MACH  card. 

42-44  IA(13) 

45-47  IA(14) 

48-50  IA(15) 

51-53  IA(16) 

54-56  IA(17) 

57-59  IA(18) 

60-62  IA(19) 

63-65  IA(20) 

3.7  GROUP  CARD 

the  plots  produced  by  this  module  are  divided  into  seven  distinct  plot 
groups  as  described  in  Section  2.2.  Each  group  contains  related  data;  such 
as  static  stability  plots,  downwash  plots,  etc.  The  specific  plots  in  each 
plot  group  are  shown  in  Figures  4  through  10.  Plot  groups  one  through  four 


are  plotted  by  default  and  groups  five  through  eeven  at  the  users  option. 
The  GR0UP  card  performs  two  functions:  (1)  it  activates,  or  deactivates,  a 
plot  group;  and  (2)  defines  the  group  number  for  subsequent  definition  of  the 
plots  that  will  be  created  and  (optionally)  their  scaling  and  axis  labels. 
Thus,  even  if  the  group  is  nominally  active  the  group  card  must  be  used  if 
these  options  are  exercised.  The  format  of  the  GR0UP  card  and  definition 
of  its  input  data  are  given  below. 

Explanation 

Control  word,  the  five  letters  GR0UP. 

Plot  group  number,  right  justified  integer. 

NGP  ”  group  number 
If  NGP  <  -1,  the  group  is  deactivated 
If  NGP  >  1,  the  group  is  qctive  and  the 
PL0T,  SCALE,  and  LABEL 
control  cards  can  refer  to 
the  group. 

Default:  Groups  1,  2,  3,  and  4  are  active, 
groups  5,  6,  and  7  are  inactive. 

Note:  Plots  are  produced  only  from  active 
plot  groups. 

3.8  PL0T  CARD 

The  PL0T  card  must  follow  a  GR0UP  card  that  defines  an  active  plot 
group  (i.e.  1  _<  NGP  <  7).  Other  control  cards  may  be  placed  between  the 
GR0UP  and  PL0T  cards;  however,  each  PL0T  card  refers  to  the  immediately  pre¬ 
ceding  GR0UP  card.  A  PL0T  card  performs  two  functions:  (1)  it  activates, 
or  deactivates,  a  plot;  and  (2)  defines  a  plot  for  subsequent  optional  defi¬ 
nition  of  its  scales  and  axis  labels.  A  PL0T  card  must  be  supplied  for  each 
plot  for  which  SCALE  or  LABEL  control  cards  are  used.  A  definition  of  the 
nominally  active  and  inactive  plots  are  given  in  Figures  4  through  10.  For  a 
plot  to  be  drawn,  it  must  be  in  an  active  plot  group.  The  format,  of  the  PL0T 
card  and  definition  of  its  input  data  are  given  below. 


Colum" 

Code 

Explanation 

1-4 

PL0T 

Control  word,  the 

four  letters  PL0T. 

6-8 

NPLT 

Plot  number,  right 

justified  integer 

NPLT  ■  plot  number  in  group  NGP  (see 
GROUP  card). 


Column  Code 

1-5  GR0UP 

7-8  NGP 
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If  NPLT  <  -1,  the  plot  is  deactivated 
If  NPLT,  2  1 »  the  plot  is  active  and  the 
SCALE  and  LABEL  control 
cards  can  refer  to  it. 

Default:  see  Figures  4  to  10 

3.9  SCALE  CARD 

The  SCALE  card  allows  the  user  to  specify  the  upper  and  lower  limits 
of  the  plot  scales.  These  limits  may  be  adjusted  so  that  the  Increment 
between  annotations  on  the  axis  is  1,  2,  4,  5,  or  8  times  •»  power  of  ten  (see 
Section  6.9).  Any  data  that  lies  outside  of  the  final  upper  and  lower  scale 
limits  will  not  be  plotted.  If  the  axis  scales  are  not  specified,  or  set 
equal  to  "UNUSED,”  the  program  automatically  defines  the  plot  scales.  The 
SCALE  card  must  follow  a  PL0T  card  that  defines  an  active  plot.  Other 
control  cards  may  be  placed  between  the  PL0T  and  SCALE  cards,  however,  the 
SCALE  card  refers  to  the  immediately  preceding  PL0T  card.  The  format  of  the 
SCALE  card  and  definition  of  its  input  data  are  given  below. 


Column 

Code 

Explanation 

1-5 

SCALE 

Control  word,  the  five  letters  SCALE. 

7 

AXIS 

Axis  identification,  either  the  letter 

X  or  Y. 

li-20 

SL 

Scale  lower  limit,  real  number. 

21-30 

SU 

Scale  upper  limit,  real  number. 

3.10  LABEL  CARD 

The  LABEL  card  allows  the  user  to  specify  the  label  that  will  appear 
along  the  axes  of  each  plot.  The  default  axis  labels  for  each  plot  are  given 
in  Figures  4  through  10.  The  LABEL  card  must  follow  a  PL0T  card  that  defines 
an  active  plot.  Other  control  cards  may  be  placed  between  the  PL0T  and  LABEL 
cards,  however,  the  LABEL  card  refers  to  the  immediately  preceding  PL0T  card. 
The  format  of  the  LABEL  card  and  its  input  data  are  defined  below. 

Column  Code  Explanation 

1-5  LABEL  Control  word,  the  five  letters  LABEL 

7  AXIS  Axis  identification.  Input  either  the  letter 

X  or  Y. 
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Column  Code  Explanation 

9-10  NC  The  number  of  characters  in  the  axis  label, 

right  justified  integer.  If  NC  *  0,  a  label 
is  not  drawn. 

11-50  Title  Text  of  the  label.  The  characters  must  be 

left  justified  in  the  field. 

Default:  uee  Figures  4  through  11. 

3.11  GRID  CARD 

The  GRID  card  controls  the  drawing  of  a  one  inch  grid  on  the  plot 
field.  The  format  of  the  GRID  card  and  its  input  data  are  given  below. 
Column  Code  Explanation 

1-4  GRID  Control  word,  the  four  letters  GRID. 

6-8  IG  Grid  code,  right  justified  integer. 

If  IG  _<  -1  no  grid  is  drawn. 

If  IG  >  0  a  one  inch  grid  is  drawn. 

Default: 

IG  -  -1. 

3.12  SAVE  CARD 

When  the  SAVE  card  is  present  in  a  case,  all  of  the  control  cards 
input  for  that  case  are  preserved  for  use  in  the  following  case.  The 
use  of  this  card  minimizes  the  inputs  required  for  multiple  case  jobs. 
The  total  number  of  control  cards  in  consecutive  save  cases  cannot  exceed 
100,  this  includes  multiple  occurrences  of  the  same  control  card.  The  format 
of  the  SAVE  card  is  given  below. 

Column  Code  Explanation 
1-4  SAVE  Control  word ,  the  four  letters  SAVE. 

3.13  NEXT  CASE  CARD 

This  control  card  terminates  the  input  data  for  a  case  and  initiates 
plotting.  Any  control  cards  following  this  card  are  in  the  following 
case.  After  the  case  plotting  is  complete,  the  program  will  return  to  the 
input  routine  and  expect  inputs  for  another  case.  The  format  of  the  NEXT 
CASE  card  is  given  below.  >  . 

Column  Code  Explanation 
1-9  NEXT  CASE  Control  word ,  the  two  word  NEXT  and  CASE 
separated  by  a  blank. 
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3.14  END  CARD 


The  END  card  performs  the  same  functions  as  the  NEXT  CASE  card;  however, 
after  the  case  plotting  is  complete,  the  program  execution  is  terminated.  An 
end-of-file  on  the  input  stream  has  the  same  effect  as  the  END  card.  The 
format  of  the  END  card  is  given  below. 

Column  Code  Explanation 
1-3  END  Control  word,  the  three  letters  END. 
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ID 


IDllllNEllSYM 


-  FILE  IDENTIFICATION  (A  OR  B) 

I UNIT  -  FILE  UNIT  NUMBER  (11  to  14  FOR 
ACTIVE  FILES  ),  INTEGER 

ICASE  -  (ASE  NUMBER  ON  FILE,  INTEGER 

ILINE  -  LINE  CODE,  INTEGER 

*  -1 ,  NO  LINE 

*  0,  STRAIGHT  LINE,  DEFAULT 

*  1,  SPLINE  CURVE 

ISVM  -  SYMBOL  CODE,  INTEGER 

=  0,  NO  SYMBOLS,  DLFAULT 

*  j.  SYMBOL  ON  EVERY  jth  point. 

NC  -  NUMBER  OF  CHARACTERS  IN  TITLE,  INTEGER 

IM  -  MACH  NUMBER  INDICES,  INTEGER 

IA  -  ALTITUDE  INDICES,  INTEGER 

AXIS.  -  AXIS  IDENTIFICATION  (X  OR  Y) 

SL  -  SCALE  LOWER  LIMIT,  REAL 

SU  -  SCALE  UPPER  LIMIT,  REAL 

IG  -  GRID  CODE,  INTEGER 

<  -1  NO  GRID  ON  PLOT,  DEFAULT 
>  0  GRIO  WILL  BE  DRAWN  ON  PLOT 


LIKE 


T  I  I  I  E  Ipj 
T.  l  T  l  E  21.1 


BJJLLDl 


i  =  i  2 


jLiia 


NGP 


'GR0.ll.Pl 


LfLLQ-T  £  j 


AXIS  NC 


[label! 


lSAV.Ej 


FIGURE  3  PLOT  MODULE  CONTROL  CARD  SUMMARY 


NOTES:  Leave  Unused  Columns  Blank 

Integer  variables  are  right  justified  In  the  Input  field. 
Titles  are  left  justified  In  the  Input  field. 

Real  variables  must  have  a  decimal  point. 
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PLOT  GROUP  1 
STATIC  STABILITY  DATA 


NO. 

TYPE* 

DESCRIPTION 

AXIS 

AXIS  LABEL 

1 

Cq  VS.  cc 

Y 

DRAG  COEFFICIENT 

A 

X 

ANGIE  OF  ATTACK. DEG 

2 

A 

CLVS.« 

Y 

LIFT  COEFFICIENT 

X 

ANGLE  OF  ATTACK.  DEG 

CmVS.<* 

Y 

PITCHING  MOMENT  COEFFICIENT 

1  « 

A 

X 

ANGLE  OF  ATTACK. OEG 

i 

Cg  VS.3 

Y 

NORMAL  FORCE  COEFFICIENT 

X 

ANGLE  OF  ATTACK.  DEG 

5 

i 

CAVS.« 

Y 

AXIAL  FORCE  COEFFICIENT 

X 

ANGLE  OF  ATTACK.  DEG  . 

s 

i 

C|yOC  VS.  3 

Y 

LIFT  CURVE  SLOPE.  1/OEG 

X 

ANGLE  OF  ATTACK.  DEG 

i 

Cm«VS.« 

Y 

PITCHING  MOMENT  COEF.  SLOPE.  1/DEG 

7 

.  him  — 

X 

ANGLE  OF  ATTACK, DEG 

j| 

i 

cY^vs.« 

Y 

SIDE  FORCE  COEF.  OERiV.,  1/DEG 

X 

ANGLE  OF  ATTACK.  DEG 

Q 

C"0VS-a 

Y 

YAWING  MOMENT  COEF.  DERIV.,  1/DEG 

1 

~  X 

'  ANGLE  OF  ATTACK.  DEG 

10 

i 

Y 

ROLLING  MOMENT  COEF.  DERIV.,  1/OEG 

X 

ANGLE  OF  ATTACK.  DEG  i 

•NOMINAL  PLOTTING  TYPE 
A  ACTIVE  PLOT 
I  -  INACTIVE  PLOT 


FIGURE  4  GROUP  1  PLOTS 
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PLOT  GROUP  2 

DRAG  POLARS  AND  STABILITY  PLOTS 


NO. 

TYPE* 

DESCRIPTION 

AXIS 

AXIS  LABEL  1 

1 

Y 

LIFT  COEFFICIENT 

A 

cLvs.  Cm 

PITCHING  MOMENT  COEFFICIENT 

cLvs.c0 

Y 

LIFT  COEFFICIENT 

2 

A 

X 

ORAG  COEFFICIENT 

cNvs.cA 

Y 

NORMAL  FORCE  COEFFICIENT 

3 

1 

X 

AXIAL  FORCE  COEFFICIENT 

A 

1 

C N  VS.  Cm 

Y 

NORMAL  FORCE  COEFFICIENT 

X  ' 

PITCHING  MOMENT  COEFFICIENT 

‘NOMINAL  PLOTTING  TYPE 
A -ACTIVE  PLOT 
I  -  INACTIVE  PLOT 


FIGURE  5  GROUP  2  PLOTS 


PLOT  GROUP  3 
DOWNWASH  DATA 


NO. 

TYPE* 

DESCRIPTION 

AXIS  LABEL 

1 

D 

VS.« 

8 

DYNAMIC  PRESSURE  RATIO 

ANGLE  OF  ATTACK, DEG 

2 

n 

e  VS;« 

8 

DOWNWASH  ANGLE,  DEG 

ANGLE  OF  ATTACK,  OEG 

3 

D 

—  vs.« 

da 

H 

DOWNWASH  GRADIENT 

ANGLE  OF  ATTACK,  DEG 

D 

a 

w,vs.« 

8 

CARNARD  EFF.  DOWNWASH  ANGLE,  OEG 
ANGLE  OF  ATTACK,  DEG 

5 

a 

(de/doc)8  VS,  “ 

8 

CARNARD  EFF.  DOWNWASH  GREDIENT 
ANGLE  OF  ATTACK,  DEG 

•NOMINAL  PLOTTING  TYPE 
A -ACTIVE  PLOT 
I  -  INACTIVE  PLOT 


FIGURES  GROUP 3 PLOTS 
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PLOT  GROUP  4 


MACH  EFFECTS 


NO. 

TYPE* 

DESCRIPTION 

AXIS 

AXIS  LABEL 

1 

A 

C0flVS.M 

_  Y_ 

X 

-  i|R0  UFT£RAG  CQE_FF!CIENJ  _  __  _ 
"  MACH  NUMBER 

2 

A 

C^VS.M 

_  Y_^ 
X 

LIFT  CURVE  SLOPE.  1/DEG 
,  MACH  NUMBER 

3 

A 

Cm<xVS.M 

_Y_ 

X 

_  PJTCH^NG_MOMENT_COEF.  SLOPE,  J/D  EG 
MACH  NUMBER 

4 

1 

Cy0vs.m 

Y 

X 

SIDE  FORCE  COEF.  DERIV.,  1/DEG 

MACH  NUMBER 

5 

1 

Cn^  VS.  M 

_  Y__ 
~  X 

YAWING  MOMENT  COEF.  DERIV..  1/OEG 
MACH  NUMBER 

0 

1 

C^VS.M 

Y 

X 

JIOLUNG^MOMENTJCOEF.  DERIV.,  1/DEG 
~  MACH  NUMBER 

•NOMINAL  PLOTTING  TYPE 
A  -ACTIVE  PLOT 
I  •  INACTIVE  PLOT 


FIGURE  7  GROUP  4  PLOTS 


PLOT  GROUP  5 
SYMETRICAL  FLAP  DATA 


NO. 

TYPE* 

DESCRIPTION 

AXIS 

AXIS  LABEL 

1 

1 

AClVS.8 

Y 

X 

LIFT  COEFFICIENT  :MREMENT_  __  __  __ 
"FLAP  DEFLECTION,  DEG 

2 

1 

ACm  VS.  8 

_  Y_  ^ 
X 

PITCHING  MOMENT  COEF.  INCREMENT 

FLAP  DEFLECTION,  DEG 

3 

1 

ACi  VS.  6 

Lnnx 

_Y_ 

X 

MAX.  LIFT  COEF.  INCREMENT 

FLAP  DEFLECTION,  DEG 

4 

1 

ACn  .  VS.  5 
urmn 

Y 

X 

MIN.  DRAG  COEF.  INCREMENT 

FLAP  DEFLECTION.  DEG 

5 

1 

ic^avs.8 

Y 

X 

(CLA)D  -  DERIV.  OF  LIFT  COEF.  SLOPE 

FLAP  DEFLECTION,  DEG 

6 

1 

chaVS.« 

Y 

X 

_CHA-_HINGEMOMENT JDERIV.J/DEG_  _ 
‘FLAP  DEFLECTION,  DEG 

7 

1 

ch8VS.« 

_Y_ 

X 

CHO  -  HINGE  MOMENT  DERIV.,  1/DEG 

FLAP  DEFLECTION,  DEG 

8 

1 

ACq.  VS.  <*,  AT  CONST. 8 

Y 

X  ~ 

INDUCED  DRAG  INCREMENT 

ANGLE  OF  ATTACK,  DEG 

•NOMINAL  PLOTTING  TYPE 


A  •  ACTIVE  PLOT 
I  •  INACTIVE  PLOT 


FIGURE  8  GROUP  5  PLOTS 


PLOT  GROUP  6 
ASYMMETRICAL  FLAP  DATA 


NO. 

TYPE* 

DESCRIPTION 

AXIS 

AXIS  LABEL 

1 

Cf  VS.  5s/C 

Y 

ROLLING  MOMENT  COEFFICIENT 

1 

X 

SPOILER  HEIGHT,  FRACTION  OF  CHORD 

* 

1 

C„VS.&/C 

Y 

YAWING  MOMENT  COEFFICIENT 

X 

SPOILER  HEIGHT,  FRACTION  OF  CHORD 

CnVS.  «  AT  CONST.  !6L-SR) 

Y 

YAWING  MOMENT  COEFFICIENT 

3 

1 

X 

ANGLE  OF  ATTACK,  DEG 

1 

'  s 

CjVS.  («,,  -«„) 

Y 

ROLLING  MOMENT  COEFFICIENT 

J 

"  "x 

TncFemThTAlTlAP  GEFUL^RirDEG 

5 

1 

CfVS.«  AT  CONST.  (5l  -5r) 

Y 

ROLLING  MOMENT  COEFFICIENT 

X 

ANGLE  OF  ATTACK,  DEG 

•NOMINAL  PLOT  TYPE 
A  •  ACTIVE  PLOT 
I  -  INACTIVE  PLOT 


FIGURE  9  GROUP  6  PLOTS 
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PLOT  GR0UP7 
TRIM  DATA 


NO. 

TYPE* 

DESCRIPTION 

AXIS 

AXIS  LABEL 

1 

■ 

- 

6,  VS. « 

Y 

TRIM  FLAP  OEFLECTION,  OEG 

X 

ANGLE  OF  ATTACK,  DEG 

mm 

I 

Y 

LIFT  COEFFICIENT  INCREMENT 

4 

■■ 

£ 

X 

ANGLE  OF  ATTACK,  DEG 

1 

< 

-J 

mm 

MAX.  LIFT  COEF.  INCREMENT 

X 

BS 

ANGLE  OF  ATTACK,  DEG 

■ 

AC0.VS.« 

§t 

mm 

INDUCED  DRAG  INCREMENT 

n 

s 

SB 

ANGLE  OF  ATTACK,  DEG 

ACn  •  VS.  a 

umin 

■n 

MIN.  DRAG  COEF.  INCREMENT 

5 

mm 

ANGLE  OF  ATTACK.  DEG 

6 

ch8VS.« 

Y 

CHD  -  HINGE  MOMENT  DERIV.,  1/DEG 

X 

ang;  e  of  attack,  deg 

a 

<J£jt  VS.  « 

Y 

TRIM  HORIZ.  TAIL  INCIDENCE,  DEG 

tS-j 

oo  <2 

<»-:= 

ujz® 

>Oh 

ONJ 

X 

ANGLE  OF  ATTACK,  DEG 

a 

CL  VS.<* 

8 

_  _UFTCOEFFICIENT  _ 

"  ANGLE  OF  ATTACK,  DEG 

a 

c0  VS. a 

JOh 

mem 

DRAG  COEFFICIENT 

<s 

8 

ANGLE  OF  ATTACK,  DEG 

•NOMINAL  PLOT  TYPE 
A  ACTIVE  PLOT 
I  -  INACTIVE  PLOT 


FIGURE  10  GROUP  7  PLOT 
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SECTION  A 


SAMPLE  CASES 


The  example  problem  discussed  In  this  section  was  designed  to  illustrate 
the  plot  module  features  and  to  use  of  all  the  available  control  cards.  Two 
sets  of  Digital  Datcom  plot  data  are  utilized  and  the  input  data  for  both 
Digital  Datcom  problems  are  shown  in  Figure  11.  Problem  one  uses  the  con¬ 
figuration  of  example  problem  3  in  Volume  I  and  is  run  at  six  subsonic  Mach 
numbers.  Case  one  is  the  basic  configuration  data  and  case  two  adds  a 
symmetrical  flap  on  the  horizontal  tail.  Problem  two  uses  the  configuration 
from  problem  one  with  the  wing  span  reduced  by  16  percent.  Case  one  is  the 
basic  configuration  data  and  case  two  adds  trim  with  an  all  moveable  hori¬ 
zontal  tail.  These  problems  were  designed  to  provide  two  data  files  with  two 
cases  on  each  file.  The  plot  module  accesses  problem  1  on  unit  11  and 
problem  2  on  unit  12. 

The  plot  module  sample  problem  consists  of  six  cases.  These  cases  were 
designed  to  illustrate  the  use  of  all  the  plot  module  control  cards  and  the 
various  program  features.  The  input  data  are  shown  in  Figure  12  and  a 
description  of  each  case  is  given  below. 

A. 1  PLOT  CASE  1 


This  case  presents  the  basic  stability  data  for  Digital  Datcom  problem 
one .  It  shows  the  plots  that  are  generated  by  default ,  form  two  of  the  BUILD 
card,  and  the  MACH  card.  The  program  will  select  Digital  Datcom  case  one  on 
1/0  unit  11  by  default,  this  is  problem  one  case  one.  This  case  has  six 
Mach  numbers,  therefore  the  MACH  card  is  used  to  select  just  one  Mach  number 
when  plotting  data  from  plot  groups,  1,  2,  and  3.  All  six  of  the  Mach 
numbers  are  used  when  plotting  the  Mach  effects,  plot  group  A.  The  11  plots 
for  this  case  are  shown  in  Figure  13. 

A. 2  PLOT  CASE  2 


This  case  presents  a  comparison  of  data  from  Vigital  Datcom  problems-  one 
and  two.  Problem  one  is  assigned  to  file  A,  I/O  unit  11,  and  problem  two  to 
file  B,  I/O  unit  12.  This  case  demonstrates  two  line  codes,  the  deactivation 
of  plot  groups,  form  one  of  the  BUILD  card,  and  use  of  the  SAVE  card.  The 
MACH  and  TITLE  cards  are  repeated  for  this  case  because  they  were  not  saved 
in  case  one.  Only  the  three  nominal  plots  from  plot  group  one  are  produced 
for  this  case  and  are  presented  in  Figure  1A. 

31 


4.3  PLOT  CASE  3 

This  case  presents  a  drag  comparison  at  low  angles  of  attack.  The  input 
data  from  case  two  were  saved,  so  only  the  additional  inputs  to  select  the 
drag  coefficient  plot  and  to  scale  the  data  are  required.  A  grid  is  speci¬ 
fied  and  the  label  of  abscissa  (Y-axis)  is  modified.  The  plot  for  this  case 
is  presented  in  Figure  15. 

4.4  PLOT  CASE  4 

This  case  presents  a  complete  configuration  build-up  cf  the  zero  lift 
drag  as  a  function  of  Mach  number.  The  input  data  fro.:  case  three  were  saved, 
which  includes  the  case  two  data.  Plot  one  cf  group  four  is  the  only  active 
plot  and  use  of  form  four  of  the  BUILD  card  is  demonstrated.  The  plot  for 
this  case  is  presented  in  Figure  16. 

4.5  PLOT  CASE  5 

This  case  presents  the  basic  stability  (group  1)  and  all  moveable 
horizontal  tail  trim  data  found  in  Digital  Datcom  problem  two,  case  two. 
The  input  data  from  case  four  was  saved  so  file  B  must  be  deactivated  and 
file  A  defined  on  1/0  unit  12,  case  two.  Plot  group  seven  is  activated  and 
plots  7,8,  and  9  are  activated  since  all  plots  in  this  group  are  nominally 
inactive.  Plot  group  one  is  reactivated  and  plots  two  and  three  are  also 
reactivated.  In  case  three,  the  scale  limits  and  abscissa  label  were  defined 
for  plot  one  of  group  one.  In  this  case  the  axis  label  is  restored  to  its 
default  value  and  automatic  scaling  is  restored  by  inserting  the  program 
vulue  for  "UNUSED"  on  the  scale  cards.  The  plot  grid  is  removed,  form  three 
of  the  BUILD  card,  and  a  spline  curve  through  the  data  points  are  demonstrat¬ 
ed.  The  six  plots  for  this  case  are  presented  in  Figure  17. 

4.6  PLOT  CASE  6 

This  case  presents  the  incremental  effects  on  lift,  drag  and  pitching 
moment  due  to  a  plain  trailing  edge  flap.  Digital  Datcom  problem  one,  case 
two.  A  complete  set  oi  input  data  are  required  because  a  SAVE  card  did  not, 
appear  in  case  five.  Plots  1  and  2  in  group  five  are  active  in  this  case, 
symbols  are  placed  on  every  second  data  point,  and  a  grid  is  drawn  on 
the  plots.  Most  computers  interpret  a  blank  input  field  as  a  zero,  therefore 
the  zero  on  the  GRID  card  is  not  required.  The  two  plots  for  this  case  are 
presented  in  Figure  18. 
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CASEID  PLOT  TEST  PROBLEM  2,  CASE  2 
PLOT 

NEXT  CASE  . 


T 


PROBLEM  1 


PROBLEM  2 


FIGURE  11  DIGITAL  DATCOM  PLOT  PROBLEMS 
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TITLEl  26  PLOT  MODULE  SAMPLE  PRC3LEM" 

TITLE2  6  CASE  I 

MACH  5 

BUILD  ALL 

NEXT  CASE 

TITLEl  26  PLOT  MODULE  SAMPLE  PRGBLEM“ 

TITLE2  6  CASE  2 

FILE  B  12  1 

LINE  A  1  0 

LINE  B  -1  1 

GROUP  -2 

GROUP  -3 

GROUP  -4 

MACH  5 

BUILD  W,  BWHV 

SAVE 

NEXT  CASE 
TITLE 2  6  CASE  3" 

GROUP  1 

PT  r>v  _2 

PLOT  -3 
PLOT  1 

SCALE  X  -2.0  12.0 

SCALE  Y  0.0  0.12 

LABEL  Y  2CD 
GRID  0 
SAVE 

NEXT  CASE  _ _ 

TITLE2  6  CASE  4 

BUILD 

GROUP  -1 

GROUP  4 

PLOT  -2 

PLOT  -3 

SAVE 

NEXT  CASE 


CASE  1 


CASE  2 


CASE  3 


CASE  4 


F1CURE  12  SAMPLE  PLOT  INPUT  DATA 


J 


IFT  COEFFICIENT  ORRC  COEFFICIENT 

.KO  .to  1.20  1.60  .  00  .OS  .10  .IS 


TCHiNC  MOMENT  COEFFICIENT 


PLOT  MODULE  SAMPLE  PROBLEM 
CASE  1  • 

MACH  NO.  =  .50  REYNOLDS  NO.  =  3.50  *10® 


=  2.250  FT  CBAR  =  .822  FT  SPAN  =  3.000  FT 
XCG  =  2.60  FT  2CG  =  .00  FT 


PLOT  MOOULE  SAMPLE  PROBLEM 
CASE  1 

MACH  NO.  =  .50  REYNOLDS  NO.  =  3.50  -10® 

=  2.250  FT.  CBAR  «  .822  FT  SPAN  =  3.000  FT 
XCG  =  2.60  FT  ZCG  =  .00  FT 


FIGURE  13  SAMPLE  CASE  1  PLOTS  (CONTINUED) 


ZERO  LIFT  ORAG  COEFFICIENT  »10 


PLOT  MODULE  SAMPLE  PROBLEM 
CASE  i 

ALTITUDE  =0  FT 

F  =  2.250  FT  CBAR  =  .822  FT  SPAN  =  3.000 
o  XCG  =  2.60  FT  ZCG  =  .00  FT 


FT 


8-U 

.oo 


.08 


.16 


.24 


.32 

MPCH  NUMBER 


.10 


.48 


.56 


H  BWHV 

y  bwv 

Z  BWH 
X  BV 

♦  BH 

❖  BW 
X  VT 
+  HT 

A  RING 
O  B0DT 

-i 

.64 


PLOT  MODULE  SAMPLE  PROBLEM 
CASE  1 

ALTITUDE  =0  FT 

=  2.250  FT  CBAR  =  .822  FT  SPAN  =  3.000  FT 


FIGURE  13  SAMPLE  CASE  1  PLOTS  (CONTINUED) 
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LIFT  COEFFICIENT  DRAG  COEFFICIENT 


PLOT  MODULE  SAMPLE  PROBLEM 
CASE  2 

MfiCH  NO.  =  .50  REYNOLDS  NO.  =  3.50  *106 

=  2.250  FT  C8flR  =  .822  FT  SPAN  =  3.000  FT 
XCG  =  2.60  FT  ZCG  =  .00  FT 


a  bwhv 
Q  WING 


8.00  12.00  16.00 
ANGLE  OF  ATTACK.  DEG 


PLOT  MODULE  SAMPLE  PROBLEM 
CASE  2 

MPCH  NO.  =  .50  REYNOLDS  NO.  =  3.50  *106 

:  2.250  FT  CBflR  =  .822  FT  SPAN  »  3.000  FT 
XCG  =  2. 60  FT  ZCG  =  .00  FT 


A  BWHV 
©  WING 


'1.00  8.00  12.00  16.00  20.00 
ANGLE  OF  ATTACK.  DEG 


2K.00  28.00 


FIGURE  14  SAMPLE  CASE  2  PLOTS 
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PLOT  MODULE  SAMPLE  PROBLEM 
CASE  3 

MACH  NO.  =  .50  REYNOLDS  NO.  *  3.50  **10® 

SREF  =  2.250  FT  ,CBAR  =  .822  FT  SPAN  =  3.000  FT 
XCG  =  2.60  FT  ZCG  =  .00  FT 


FIGURE  15  SAMPLE  CASE  3  PLOT 
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ZERO  LIFT  DRAG  COEFFICIENT  -lO'1 


PLOT  MOOULE  SAMPLE  PROBLEM 
CASE  4 

ALTITUDE  =0  FT 

SREF  =  2.250  FT  CBAR  =  .822  FT  SPAN  =  3.000  FT 
XCG  =  2.60  FT  ZCG  =  .00  FT 
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X  BWH 

+  BN 
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O  BOOT 
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1 3. 

A 

A 

A 

A 

A 
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*  1 

00 

08 

16 

32 

40  . 

48  .! 

56  .1 

MPCH  NUMBER 


FIGURE  16  SAMPLE  CASE  4  PLOT 


PLOT  MODULE  SAMPLE  PROBLEM 
CASE  5 

MflCH  NO.  =  .50  REYNOLDS  NO.  =  3.50  *10® 

SREF  =  2.250  FT  CBRR  =  .822  FT  SPAN  =  3.000  FT 
XCG  =  2.60  FT  ZCC  =  .  HO  FT 


CIO  .00  4.C3  8.00  12.00  16.00  20.00  24.00 

ANGLE  OF  ATTACK.  OEG 

PLOT  MODULE  SAMPLE  PROBLEM 
CASE  5 

MflCH  NO.  *  .50  REYNOLOS  NO.  =  3.50  *10® 

SREF  =  2.250  FT  .  CBflR  =  .822  FT  SPRN  =  3.000  FT 
XCG  =  2.60  FT  ZCG  =  .00  FT 


8.00  12.00  16.00  20.00 
ANGLE  OF  ATTACK.  DEG 


FIGURE  17  SAMPLE  CASE  5  PLOTS  (CONTINUED) 


PLOT  MODULE  SRMPLE  PROBLEM 
CASE  6 

MACH  NO.  =  .50  REYNOLDS  NO.  =  3.50  «106 

SREF  =  2.250  FT  C8AR  -  .822  FT  SPAN  =  3.000  FT 
XCG  =  2.60  FT  ZCG  =  .00  FT 


4.00  8.00 


12.00  IS. 00  20.00  2H.Q0  28.00  32.00 

FLOP  DEFLECTION.  DEG 


PLOT  MODULE  SRMPLE  PROBLEM 
CASE  6 
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SECTION  5 


PROGRAM  STRUCTURE 

The  Digital  Datcom  plot  module  was  written  in  Fortran  IV  for  the  CDC 
Cyber  74  and  175  computer  systems.  Users  should  refer  to  Section  6  of  Volume 
II  when  implementing  this  module  on  another  computer  system. 

5.1  SYSTEM  RESOURCE  REQUIREMENTS 

The  plot  module  requires  the  availability  of  certain  computer  resources 
and  capabilities  in  order  to  operate.  These  resources  are: 

o  Seven  disk  or  tape  files  are  required  for  manipulation  and  processing 
of  data.  The  logical  I/O  units  are  8  through  14  and  are  in  addi¬ 
tion  to  the  system  input  and  output,  assigned  to  I/O  units  5  and  6 
respectively. 

o  The  system  should  have  an  overlay  capability,  however,  the  plot 
module  could  operate  without  this  capability, 
o  The  system  must  have  a  Fortran  compiler  that  provides  statement 
transfer  when  an  end  of  file  is  encountered. 

Each  logical  1/4  unit  used  by  this  module  is  reserved  for  a  specific 
purpose.  These  units  and  their  use  by  the  program  are  listed  below. 

Unit  Program  Usage 

5  Standard  system  input  (card  reader). 

6  Standard  system  output  (printer). 

8  Temporary  stprage  for  input  data. 

9  Temporary  storage  for  Mach  effects  data. 

10  Plot  tape. 

11  Digital  Datcom  plot  data  file. 

12  Digital  Datcom  plot  data  file. 

13  Digital  Datcom  plot  data  fiie. 

14  Digital  Datcom  plot  data  file. 

5.2  PROCRAM  DESCRIPTION 

The  plot  module  overlay  structure  Is  composed  of  a  root  segment  (overlay 
0)  and  ten  primary  overlays.  Table  1  shows  the  overall  program  structure  and 


lists  the  routines  in  each  overlay.  Table  2  lists  each  routine,  the  overlay 
where  it  is  referenced,  and  a  brief  description  of  its  purpose.  A  complete 
listing  of  each  routine  is  provided  as  a  microfiche  supplement  to  this 
report. 

5.3  DIGITAL  DATCOM  PLOT  FILE  STRUCTURE 

A  Digital  Datcom  plot  file  is  written  to  1/0  unit  13  when  the  PL0T 
control  card  is  present  in  a  case,  see  Volume  I.  This  file  contains  Digital 
Datcom  case  data  (described  below)  and  aerodynamic  data  for  each  Mach  number 
and  altitude  combination  computed  by  the  program.  The  aerodynamic  data  con¬ 
sists  of  the  first  200  elements  of  each  Ideal  Output  Matrix  (IOM)  data  block, 
the  zero  lift  drag  of  each  component  combination  in  the  IOM,  and  flap  input 
and  IOM  arrays  when  flap  or  trim  data  are  calculated.  If  the  user  generates 
a  plot  file  containing  data  from  wind  tunnel  test  of  other  sources,  those 
elements  in  the  file  that  are  defined  should  contain  the  program  value  of 
UNUSED  (1  x  10  ^  on  CDC  systems).  The  specific  data  in  the  Digital  Datcom 
plot  file  and  the  format  used  to  write  it  are  given  below. 

5.3.1  Digital  Datcom  Case  Data 

The  Digital  Datcom  case  data  consists  of  the  flight  condition  inputs 
(160  elements),  synthesis  inputs  (15  elements),  reference  parameter  inputs 
(A  elements),  and  the  logical  flags  in  the  Digital  Datcom  common  block.  FL0L0G 
(34  elements).  A  definitioin  of  the  data  in  common  block  FL0L0C  is  given  in 
Volume  II,  Section  7,  Table  10.  The  flight  condition  inputs  are  supplied  to 
Digital  Datcom  through  namelist  FLTC0N  and  are  defined  in  Volume  I,  Figure  3. 
The  synthesis  parameters  are  input  through  namelist  SYNTHS  (Volume  I,  Figure 
5)  and  the  reference  parameters  through  namelist  0PTINS  (Volume  I,  Figure  4). 
The  common  block  FL0L0C  is  defined  in  Volume  II,  Section  7.  All  of  the  input 
data  are  written  with  the  format  “1P10E12.4"  and  the  logical  flags  (FL0L0G) 
are  written  with  the  format  “34L3".  The  case  data  then  consists  of: 

1.  flight  condition  data,  16  lines  with  10  elements  on  each  line; 

2.  synthesis  data,  2  lines  with  10  elements  on  the  first  line  and  5  on 
the  second; 

3«  reference  parameters,  1  line  with  4  elements;  and 

4.  Digital  Datcom  logical  flags,  1  line  with  34  elements. 

5.3.2  Digital  Datcom  Aerodynamic  Data 

One  set  of  aerodynamic  data  is  written  for  every  Mach  number  -  altitude 
combination  executed  by  Digital  Datcom.  Each  set  of  aercUynamic  data 
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consists  of  the  basic  stability  data,  zero  lift  drag  data,  and  flap  data  if 
flap  or  trim  computations  are  performed.  All  of  these  data  are  written  using 
the  format  “ IP  10E12 .4”.  The  basic  stability  data  consists  of  the  first  200 
elements  of  the  11  IOM  common  blocks  IB0DY ,  IWING,  IHT,  1VT,  IVF,  IBW,  IBH, 
IBV,  IBWll,  IBWV,  and  IBWHV  plus  100  elements  of  downvash  data  from  10M  common 
block  1DWASH  and  the  canard  effective  t’ownwash  data  in  common  block  FACT. 
Each  of  these  common  blocks  are  defined  in  Volume  I,  Appendix  C.  The  zero 
lift  drag  data  consists  of  one  value  for  each  of  the  11  component  combina¬ 
tions  in  the  IOM.  The  flap  data  consists  of  the  final  200  elements  of  the 
IOM  common  blocks  IB0DY  and  IWING  plus  the  final  180  elements  of  the  blocks 
IHT  and  IVT.  One  set  of  aerodynamic  data  will  consist  of: 

1.  basic  stability  data,  230  lines  with  10  elements  on  each  line; 

2.  zero  lift  drag  data,  2  lines  with  10  elements  on  the  first  line; 
and 

3.  flap  data  if  trim,  symmetrical  flaps  or  asymmetrical  flaps  are 
present,  76  lines  with  10  elements  on  each  line. 


TABLE  1  DIGITAL  DATCOM  PLOT  MODULE  OVERLAY  DESCRIPTION 


TABLE  1  DIGITAL  DATCOM  PLOT  MODULE  OVERLAY  DESCRIPTION 


TABLE  2  DIGITAL  DATCOM  PLOT  MODULE  ROUTINE  DESCRIPTION 
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APPENDIX 


ASD  CALCOMP  LIBRARY  CALLS 


Each  computer  center  with  graphics  capability  has  a  unique  software 
package.  No  industry  standard  is  available,  however,  the  CALCOMP  library,  or 
a  derivative,  is  available  at  most  computer  centers.  The  Digital  Datcom  Plot 
Module  was  written  to  utilize  the  WP-AFB  Aeronautical  Systems  Division  (ASD) 
Computer  Center  CALCOMP  software  library.  In  order  for  the  plot  mtdule  to 
be  adapted  to  another  software  package,  each  of  the  ASD  CALCOMP  library  sub¬ 
routines  called  by  the  program,  as  defined  in  this  section,  must  be  replaced 
by  the  corresponding  references  to  the  new  library  package.  These  defini¬ 
tions  were  taken  from  the  ASD  Computer  Center  CALCOMP  Plotter  Guide.  Along 
with  each  subroutine  description,  the  Plot  Module  subroutines  from  which  it 
is  called  and  the  number  of  calls  in  that  subroutine  are  given. 

1.  AXIS  -  AXIS  DRAWING  ROUTINE 


Purpose :  This  routine  will  draw  an  axis  with  tic.  marks.  Each  tic  mark 

is  annotated  and  an  axis  title  is  drawn. 


Calling  Sequence: 


CALL  AXIS  (X,Y, BCD, NC, SIZE, THETA,  AMIN, DA) 

where  (X,Y)  The  coordinates  of  starting  point  of  the  axis,  in 
inches,  from  the  plot  reference  point  (or  origin). 

BCD  -  An  array  containing  the  title  of  the  axis.  This 
title  will  be  centered  along  the  axis  and  drawn 
parallel  to  it. 


NC 


-  The  number  of  characters  in  the  title.  If  NC 
is  positive,  annotation  will  be  on  the  counter¬ 
clockwise  side  of  the  axis.  If  NC  is  negative, 
annotation  will  be  on  the  clockwise  side  of  the 
axis. 


SIZE  -  The  length  of  the  axis  to  be  drawn,  In  Inches. 

THETA  -  The  anp.le  of  the  axis  to  be  drawn  (i.e.,  0.0  for  the 
X-axis  and  90.0  for  the  Y-axis). 

AMIN  -  The  adjusted  minimum  of  a  data  array.  When  subrou¬ 
tine  SCALE  is  used,  AMIN  should  be  A(N*K+1).  See 
Section  9  of  the  Appendix  for  a  definition  of  sub¬ 
routine  SCALE  and  the  array  A. 

DA  -  The  adjusted  value  or  length  of  the  variable  per 
inch.  When  subroutine  SCALE  is  used,  DA  should  be 
A(N*K+K+1). 

Comments:  Using  the  array  A  in  the  example  used  for  subroutine  SCALE, 

the  following  call  will  generate  an  axis  in  the  X-direction 
with  the  proper  scale. 

CALL  AXIS  (0.,0.,6HX-AXIS,-6,8.,0.,A(81),A(82)). 

To  generate  an  axis  in  the  Y-dlrection  use 

CALL  AXIS  (0. ,0. ,6HY-AXIS,6,8. ,90. ,A(81 ) , A(82 ) ) . 

The  axis  will  start  at  the  origin  (0,0)  and  extend  8  -  inches 
in  either  the  X  or  Y  directions  with  annotated  tic  marks  and 
labeling. 

Plot  Module  References: 

Called  from  subroutine  DRAW  in  two  places. 


2.  FLINE  -  CURVE  SMOOTHING  SUBPROGRAM 


Purpose :  FLINE  draws  a  smooth  curve  through  an  array  of  points. 


ai££t'juiscisi^3 


Calling  Sequence: 


CALL  FLINE  (X,Y,N,K,  ,<,L) 

where  X  -  The  array  of  unsealed  abscissa  values.  ' 

Y  -  The  array  of  unsealed  ordinate  values. 

N  -  The  number  of  points  in  the  array. 

If  N  is  positive  this  routine  will  function  the  same  as 
subroutine  LINE. 

If  N  is  negative  the  points  will  be  connected  with 
a  smooth  curve  rather  than  a  straight  line.  A  modified 
spline-fitting  technique  is  used. 

K  -  The  repeat  cycle  of  a  mixed  array  (i.e.,  the  data  to  be 
plotted  is  stored  in  every  Kth  cell  see  Section  11  of  the 
Appendix  for  a  definition  of  the  K  parameter). 

J  -  The  plot  mode  designator: 

J  <  0,  Draw  a  symbol  at  every  Jth  point. 

J  ■  0 ,  Connect  all  points  with  a  straight  line  or 
a  smooth  curve. 

J  >  0,  Connect  all  points  with  a  straight  line  or  a 
smooth  curve  and  draw  a  symbol  at  every  Jth 
point. 

L  -  An  integer  specifying  the  symbol  to  be  used  (see  Table  3 
for  for  the  symbols).  The  integers  (0-13)  inclusive  are 
the  only  ones  that  will  produce  a  centered  symbol. 
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Comments :  The  X  and  Y  data  arrays  must  be  dimensioned  at  least 

( N*K+K+ 1 )  and  the  adjusted  XMIN,  DX,  YMIN  and  DY  must 
be  stored  in  X(N*K+1),  X(N*K+K+1),  Y(N*K+1),  and  Y(N*K+K+1) 
respectively. 

Plot  Module  References: 

Called  from  subroutine  DRAW  in  one  place. 

GRID  -  GRID  DRAWING  SUBPROGRAM 


DY  -  The  spacing  or  distance,  in  inches,  between  Y 

grid  lines. 

M  -  The  number  of  X  grid  lines  to  be  drawn. 

N  -  The  number  of  Y  grid  lines  to  be  drawn. 

Comments: 

Plot  Module  References: 

Called  from  subroutine  DRAW  in  one  place. 
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4.  ,  LINE  -  ARRAY  PLOTTING  ROUTINE. 

Purpose:  This  routine  is  used  to  plot  an  array  of  points. 

Calling  Sequence: 

CALL  LINE  (ABSC,ORD.N,K,J,L) 

where  ABSC  -  The  array  of  unsealed  abscissa  (X)  values. 

ORD  -  The  array  of  unsealed  ordinate  (Y)  values. 

N  -  The  number  of  points  to  be  plotted  from  the  array. 

K  -  The  increment  used  to  select  the  points  to  be  plotted 

from  the  array  (see  Appendix  C,  for  a  discussion  of 
the  use  of  the  K  parameter).  Normally  K  *  1. 

J  -  The  plot  mode  designator: 

J  <  0,  Draw  a  symbol  at  every  Jth  point. 

J  *  0,  Connect  all  points  with  a  straight  line. 

J  >  0,  Connect  all  points  with  a  straight  line 

and  draw  a  symbol  at  every  Jth  point. 

L  -  An  integer  specifying  the  symbol  to  be  used  (see 
Symbol  Table  in  Section  6.12). 

Comments :  The  X  and  Y  data  arrays  must  be  dimensioned  at  least 

(N*K+K+1),  and  the  adjusted  XMIN,DX,YMIN  and  DY  must  be 
stored  in  X(N*K+1),  X(N*K+K+1),  Y'N*KK)  and  Y(N*K+K+1) 
respectively. 

CALL  LINE  (X,Y, 80, 1,5,0) 


>■*4 


will  connect  the  60  points  defined  by  the  (X,Y)  coordinate 
pairs  with  straight  lines  and  plot  the  symbol  every  5th 
point.  The  scale  factors  must  be  located  in  X(81),X(82), 
Y(81 )  ind  Y(82). 

Plot  Module  References: 

Called  from  subroutine  DRAW  in  one  place. 

NUMBER  -  A  ROUTINE  TO  DRAW  FLOATING  POINT  NUMBERS. 

Purpose :  This  routine  is  used  to  annotate  a  plot  with  floating  point 

numbers;  it  will  convert  a  floating  point  number  to  its 
display  code  equivalent  and  draw  the  number. 

Calling  Sequence: 

CALL  NUMBER  (X,Y,HGHT,FPN, THETA, N) 

where  (X,Y)  -  The  page  coordinates  of  the  lower  left  corner 
of  the  first  annotation  character. 

HGHT  -  The  height,  in  inches,  of  the  number  to  be  drawn. 

FPN  -  The  number  to  be  drawn. 

THETA  -  The  angle  of  lettering  with  respect  r.o  the  X-axis. 

N  -  Decimal  character  count 

N  >  0,  N  is  the  number  of  plotted  digits  to 
the  right  of  the  decimal  point. 

N  *  0,  No  fractional  part. 

N  .»  -1,  No  fractional  part  and  no  decimal  point. 


64 


N  <  -1,  N  -1  digits  are  truncated  starting  at  the 
units  digit. 

In  all  cases,  the  rightmost  digit  drawn 
will  be  rounded. 

Comments :  Subroutine  NUMBER  calls  subroutine  SYMBOL  to  draw  the  actual 

annotation  characters.  All  size  and  spacing  conventions  are 
the  same  as  SYMBOL. 

CALL  NUMBER  ( 5. ,5. ,0. 21 ,53. 2326,0. ,3) 

The  preceding  call  will  draw  the  number  53.233  parallel  to 
the  X-axis  with  the  lower  left  hand  corner  of  the  number  5 
starting  at  the  page  coordinates  (5., 5.).  Each  letter  will 
be  0.21  Inches  in  height  with  the  entire  number  taking  1.26 
inches. 

Plot  Module  References; 

Called  from  subroutine  DRAW  in  nine  places. 

Called  from  subroutine  DRAWS  in  one  place. 

Called  from  subroutine  DRAW6  in  one  place. 

6.  PLOT  -  PEN  MOVEMENT  ROUTINE. 

Purpose :  This  routine  causes  the  pen  to  move  from  its  current  position 

to  a  new  location,  which  may  represent  a  new  origin. 

Calling  Sequence: 

CALL  PLOT  (X.Y.IO 

where  (X,Y)  -  The  page  coordiantes,  in  inches,  of  the  point 
to  which  the  pen  is  to  be  moved. 
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Comments : 


IC  -  Pen  positioning  parameter 

(1)  The  magnitude  of  IC=M+N 

where  H  =  2  for  pen  down 

M  “  3  for  pen  uf> 

N  »  0  if  X  and  Y  are  physical  page 
coordinates  and  require  no 
scaling. 

N  =  10  if  the  scale  factors 

provided  through  subroutine 

OFFSET  are  to  be  used  to 

compute  the  pen  movement 

(i.e.,  convert  to  physical  page 

coordinates) 

N  ■  20  if  the  buffers  are  not  to  be 
cleared. 

(2)  If  IC  is  negative,  a  new  plot  reference  point 
is  established  (i.e.,  after  moving  to  (X,Y) 
in  linkage,  assign  this  point  as  the  new 
origin).  On  an  off-line  system  this  will  also 
result  in  writing  a  new  block  address  record. 

This  routine  is  the  workhorse  of  the  plotting  package; 

all  incremental  plotting  commands  are  generated  here. 

CALL  PL0T(5.,  6.,  3)  Will  move  the  pen  from  its  current 

position  to  the  point  (5.,  6.)  with  the 
pen  up. 

CALL  PL0T(5.t  6.,  -2)  Will  move  the  pen  from  its  current 

position  to  the  point  (5.  ,6.)  with  the 
pen  down  and  will  assign  this  point  as 
the  new  origin. 
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Plot  Module  References: 


Called  from  subroutine  DRAW  in  seven  places. 

Called  from  subroutine  DRAW5  in  two  places. 

Called  from  subroutine  DRAW6  in  two  places. 

Called  from  subroutine  WAIT  in  one  place. 

7..  PLOTE  -  PLOT  TERMINATION  ROUTINE. 

Purpose:  This  routine  creates  the  necessary  instructions  to  properly 

terminate  plotting.  It  will  print  a  time  estimate  for  the 
plot  (off-line  only),  create  terminal  instruction  for  the 
1  operator,  and  ascertain  that  the  last  buffer  or  partial  buffer 
has  been  written. 

Calling  Sequence: 

CALL  PLOTE ( N ) 

Comments :  This  routine  must  be  called  after  all  plotting  has  been 

completed.  It  should  appear  at  the  termination  of  the  main 
program,  immediately  preceding  the  STOP  statement.  The 
parameter  N  will  contain  the  number  of  the  last  block  written 
on  the  plot  tape. 

Plot  Module  References: 

Called  from  subroutine  INPUT  in  one  place. 

8.  PLOTS  -  PLOT  INITIALIZATION  ROUTINE. 

Purpose :  This  routine  (off-line  only)  is  used  to  initialize  the 

'  CALCOMP  plotting  package.  The  routine  allocates  buffers  and 
initializes  all  I/O. 


Calling  Sequence; 

CALL  PLOTS  (DATA.N.M) 

where  DATA  -  The  starting  address  of  an  area  of  core  storage 
reserved  as  a  work  region  for  the  CALCOMP  plotting 
package  (this  array  name  must  be  dimensioned  at 
least  N. 

N  -  The  number  of  words  in  this  work  region  (1024  is 
optimum). 

M  -  The  FORTRAN  logical  tape  number. 

Comments :  This  routine  must  be  called  before  any  other  routine  if  the 

off-line  plotting  package  is  referenced;  it  should  be  the 
first  executable  statement  of  the  main  program. 

Plot  Module  References: 

Called  from  MAIN  in  one  place. 

SCALE  -  ROUTINE  TO  ESTABLISH  SCALING  VALUES. 

Purpose :  This  routine  will  find  the  scale  values  of  an  array  and 

adjust  them  to  optimize  the  plot  but  maintain  reasonable 
values  for  the  axis  annotation. 


Calling  Sequence: 

CALL  SCALE  (A,S,N,K) 

where  A  -  The  array  to  be  scanned  for  minimum  values.  An  adjusted 
minimum  will  be  stored  in  A(N*K+1)  and  an  adjusted 
delta  or  length  of  variable  per  inch  will  be  stored  in 
A(N*K+K+1). 


S  -  The  length  of  axis  (in  inches)  over  which  the  data  is  to 
be  plotted. 

N  -  The  number  of  points  to  be  scanned  in  the  array. 

K  -  The  magnitude  of  K  is  the  increment  by  which  subroutine 
SCALE  will  scan  through  the  array  (see  Section  11  of  the 
Appendix  for  a  discussion  of  the  use  of  the  K  parameter). 

Comments :  Note  that  the  data  array  (A)  must  be  dimensioned  at  at  least 
(N*K+K+1).  The  minimum  selected  will  not  necessarily  be  the 
exact  data  minimum;  likewise,  the  delta  selected  is  limited 
to  1,2, A, 5,  or  8  times  a  power  of  ten.  These  selections  are 
made  in  an  attempt  to  make  the  axis  annotation  more  legible; 
however,  the  nature  of  the  data  may  be  such  that  it  will 
result  in  an  undesirable  compression  of  the  curve.  In  this 
event,  it  should  be  noted  that  the  use  of  subroutine  SCALE  is 
optional.  If  the  programmer  knows  the  range  of  his  data,  he 
may  supply  the  minimum  and  the  delta  himself.  This  technique 
is  employed  when  plotting  a  family  of  curves,  because  the 
same  scale  factor  must  be  applied  to  all  dependent  variables. 

DIMENSION  A(82) 

•  •  • 

CALL  SCALE  (A, 8, 80,1) 

Will  scan  A(l)  through  A(80)  for  the  minimum  value  of  A. 
The  data  will  be  scale  adjusted  for  an  8-inch  axis 
with  the  adjusted  minimum  being  returned  in  A(81)  and 
adjusted  delta  in  A(82). 

Plot  Module  References; 

Called  from  subroutine  DRAW  in  four  places. 
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10.  SYMBOL  -  SYMBOL  DRAWING  ROUTINE. 


Purpose :  This  routine  is  used  to  annotate  the  plot  or  draw  alpha¬ 

numeric  characters  or  special  symbols. 

Calling  Sequence:  , 

CALL  SYMBOL  (X, Y, HGHT, IBCD, THETA, N) 

where  (X,Y)  -  The  page  coordinates  of  the  lower  left  corner 
(before  rotation)  of  the  first  character for 
alphanumeric  characters.  For  a  special  symbol 
(N  is  negative) , the  page  coordinates  define  the 
center  of  the  symbol.  Note  that  only  the  first 
14  characters  (0-13)  are  centered  when  using 
a  special  symbol.  If  the  integer  specifying 
the  special  symbol  is  greater  than  16,  the  page 
coordinates  define  the  lower  left  corner  of  the 
character  drawn. 

HGHT  -  The  desired  letter  height,  in  inches.  For  best 
results,  HGHT  should  be  a  multiple  of  seven  times 
the  plotter  increment  size  (i.e.,  7*M*0. 005*0. 035*M 
for  off-line). 

IBCD  -  The  array  containing  the  BCD  character  string  or  an 
integer  specifying  a  special  character  (see  Table  3 
for  the  symbols). 

THETA  -  The  angle  of  the  lettering  with  respect  to  the 
X-axis . 

N  -  Character  count 

N  >  0,  N  is  the  number  of  characters  to  be  drawn 
from  the  alphanumeric  array  IBCD. 
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N  »  0 ,  One  alphanumeric  character  is  drawn. 

The  display  code  representation  of  this 
character  must  be  right- justified  in  IBCD. 

N  ■  -1,  One  alphanumeric  character  is  drawn. 

The  integer  specifying  this  character 
from  the  symbol  table  is  stored  in  IBCD. 
The  pen  is  also  raised. 

N  ■  -2,  Same  as  N  ■  -1  except  that  the  pen  is 
lowered  and  a  line  is  drawn  from  the 
current  pen  location  to  the  position 
specified  by  (X,Y). 

Comments:  s 

The  following  rules  govern  character  spacing: 

(1)  The  width  of  each  letter  is  (4/7)  *  HGHT. 

(2)  The  space  from  the  beginning  of  one  letter  to  the  beginning  of  the 
following  letter  is  the  same  as  HGHT. 

(3)  If  a  centered  symbdl  (N  is  negative)  is  specified,  the  spacing  and 
width  of  the  character  will  be  the  same  as  HGHT. 

CALL  SYMBOL  (5. ,5. ,0.21 ,6HSAMPLE,0. ,6) 

The  preceding  call  will  draw  the  word  "SAMPLE"  parallel  to  th- 
X-axis  with  the  lower  left  hand  corner  of  the  letter  S  starting  at 
the  page  coordinated  (5. ,5.).  Each  letter  will  be  0.21  inches  in 
height  with  the  entire  word  taking  1.26  inches. 
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Plot  Module  References: 


Called  from  subroutine  DRAW  in  nine  places. 

Called  from  subroutine  DRAW5  In  one  place. 

Called  from  subroutine  DRAW6  in  one  place. 

11.  ' USE  OF  THE  K  PARAMETER 

The  parameter  K,  which  defines  the  data  increment  in  subroutine  ..GALE 
and  LINE,  may  be  used  as  follows: 

A.  To  pick  out  all  points  in  a  single  dimensioned  array,  in  this 
case  K»1 . 

B.  To  thin  a  single  dimensioned  array;  that  is,  plot  every  Kth  point 
only.  The,  data  array  is  composed  of  N*K  points.  If  the  original  data  array 
contained  M  points,  then  N*K*M  and  the  array  must  be  dimensioned  at  least 
N*K+K+1  or  equivalently  M+K+l. 

C.  To  select  the  most  significant  part  from  a  double  precision  array, 
in  this  case  K-2. 

D.  To  select  a  given  variable  from  a  packed  array  (K“the  number 
of  variables  in  the  array).  Consider  the  packed  array  A: 

A( 1)  -  YI(1) 

A(2)  -  Y2(  1) 

A(3)  -  Yi(2) 

A<4)  -  Y2(2) 

A(5)  -  Yl(3) 

A(6)  -  Y2(3) 

A(7)  ■  Yl(min.) 

,  A(8)  ■  Y2(mln.) 

A(9)  -  Y1 (delta) 

A(10)-  Y2  (delta) 


In  this  example  N-3  and  K-2.  It  should  be  noted  that,  when  using 
a  packed  array,  K*(N+2)  cells  are  required  to  hold  the  data  and  associated 
scale  factors.  The  following  logic  might  be  used  to  scale  this  family  of  two 
curves : 

N-3 

K-2 

NK«N*K 

NKK-NK+K 

CALL  SCALE  (A,S,NK,1) 

A(NKK+l)»A(NK+2) 

A( NKK+2 )mA( NK+2 ) 

A( NK+2  )*A(, NK+1 ) 

The  programmer  should  also  note  that,  when  plotting  a  packed  array,  both 
the  abscissa  and  ordinate  arrays  must  be  packed  using  the  same  value  of  K. 
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